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Dispersion coefficients between the alkali metal atoms (Li-Rb) and alkaline- earth metal atoms 
(Be-Sr) are evaluated using matrix elements computed from frozen core configuration interaction 
calculations. Besides dispersion coefficients with both atoms in their respective ground states, 
dispersion coefficients are also given for the case where one atom is in its ground state and the other 
atom is in a low lying excited state. 

PACS numbers: 



I. INTRODUCTION 



II. METHODOLOGY OF CALCULATION 



Recently a systematic ab-initio investigation of the 
ground state potential curves for the MgX dimers, where 
X is an alkali atom, was reported [![. The motivation 
for this work was the growing interest in the produc- 
tion of molecules from ultracold atomic gases H. Such 
molecules can be formed by photo-association , or by 
Feshbach resonance tuning [3, Q . Most focus has been on 
diatomic molecules consisting of two alkali atoms 0413 ■ 
However it has been suggested that ultracold molecules 
in 2 S states would be good s yste ms for experiments on 
controlled chemical reactions Such molecules could 
be formed from an alkali atom in its ground state and 
an alkaline-earth atom (or Yb) in its ground state. Re- 
cently, the vibrational spectra of CaLi and SrLi were in- 
vestigated for their sensitivity to the m e /M p mass ratio 
[12]. There have been previous studies of the structure 
of alkali/alkaline-earth dimers p^| - [T8| . The most recent 
theoretical investigations of these molecules motivated by 
cold atom physics include LiSr [Hj], LiYb (HSU, RbSr 
[22| and RbYb (23U25j ]. A more comprehensive investi- 
gation has been made of the LiX dimers, where X is a 
alkaline-earth atom (2(| (27j ■ 

The present article investigates the long range interac- 
tion of various combinations of alkali and alkaline-earth 
atoms. The most efficient method for determining the 
long-range interaction is by computing the dispersion co- 
efficients since this leads to the factorization of one large 
many-body calculation into two smaller many-body cal- 
culations. The dispersion coefficient calculations were 
not restricted to the respective ground states. Disper- 
sion coefficients are also given for an alkali atom in its 
ground state and the nsnp 3 P° and nsnp 1 P° excited 
states of the alkaline-earth atoms. Coefficients are given 
for the nsnd 1 - 3 D e excited states of calcium and stron- 
tium since these nsnd 1 D e states have a smaller excita- 
tion energy than the nsnp 1 P° excited states. The dis- 
persion coefficients for the alkaline-earth atoms in their 
ground states and three of the lowest excited states of 
the alkali atoms are also given. The present work gives 
a comprehensive overview of the long range interactions 
between the ground and the low lying excited states of 
the alkali/alkaline-earth dimers. 



A. Overview of van der Waals interaction 
calculation 

The long-range van der Waals interaction between two 
hetero- nuclear atoms (i.e. the two atoms are different), 
with one atom in an S-state, can be function of 

inter-nuclear separation, R written |28l - l30| 



V(R) = -— - — - — 
K ' R e R s R 1Q 



(1) 



The C n parameters are the dispersion coefficients. 

The approach used to generate the dispersion coeffi- 
cients is based on the use of oscillator strength sum rules 
[28l |29| . This reduces the calculation of the C n parame- 
ters for two spherically symmetric atoms to summations 
over the products of the absorption oscillator strengths 
(originating in the ground state) divided by an energy de- 
nominator. The sums should include contributions from 
all discrete and continuum excitations. In practice, a 
pseudo-state representation is used which gives a discrete 
representation of the continuum (30l - [3^ |. The sum over 
oscillator strengths needs to be rewritten in terms of a 
sum over the reduced matrix elements of the electric mul- 
tipole operator in cases where one (or both) of the atoms 
is in a state with L > [30|. 

The major part of any calculation involves the genera- 
tion of the lists of reduced transition matrix elements for 
the two atomic states. This involves quite lengthy calcu- 
lations to generate the excitation spectrum of the pseudo- 
state representation. It is then a relatively straightfor- 
ward calculation to process the lists of matrix elements 
and generate the dispersion coefficients [3(J HH . 



B. Structure model: The alkali atoms 

The transition arrays for the alkali atoms are essen- 
tially those which were used in calculations of the dis- 
persion interactions between these atoms and the ground 
states of hydrogen and helium [30| . 

These were computed by diagonalizing the fixed core 
Hamiltonian in a large basis of Laguerre Type Orbitals 



2 



(LTO). The core Hamiltonian is based upon a Hartree- 
Fock (HF) description of the core with a semi-empirical 
core polarization potential tuned to reproduce the ener- 
gies of the low lying spectrum. The oscillator strengths 
(and other multipole expectation values) were computed 
with o per ators that included polarization corrections 

Core excitations are included in the C n calculations. 
Oscillator strength distributions were constructed by us- 
ing independent estimates of the core polarizabilities to 
constrain the sum rules [H], [38-40]. The methodol- 
ogy of using constrained sum rules to construct pseudo- 
oscillator strength distributions has been widely used 

a 



TABLE I: The static dipole (ai) and quadrupole (02) po- 
larizabilities (in a.u.) for the ground states of the alkali 
and alkaline-earth atoms. A recent review summarizes static 
dipole polarizabilities calculations and experiments 



Atom 


CICP 


MBPT-SD 


Other 


ai 


Li 


164.21 


164.08 [43] 


164.11(3) (Hylleraas) [44] 










164.2(1.1) Expt. [45] 


ai 


Na 


162.8 


163.0 [46] 


162.6(3) (Hybrid) [46J 










162.7(8) Expt. [47] 


ai 


K 


290.0 


289.1 [46J 


290.2(8) (Hybrid) [46J 










290.8(1.4) Expt. [48J 


Oil 


Rb 


315.7 


317.4 [46] 


318.6 (Hybrid) [46J 










318.8(1.4) Expt. [48J 


OLi 


Li 


1,424 


1,424(4) [49] 


1,423.26 (Hylleraas) [50] 


a 2 


Na 


1,879 


1,885(26) [49] 




a 2 


K 


5,005 


5,000(45) [49] 










5,018 [51] 




a 2 


Rb 


6,480 


6,520(80) [49] 








CICP 


CI+MBPT 


Other 


ai 


Be 


37.69 


37.76 [52] 


37.755 (ECG) [53] 


ai 


Mg 


71.35 


71.33 [52] 


74.9(27) (Hybrid) [54, 55] 


ai 


Ca 


159.4 


159.0 [52] 


157.1(13) (Hybrid) [52J 


ai 


Sr 


197.6 


202.0 [52] 


197.2(2) (Hybrid) [52J 








198.9 [56] 


197.14(20) (Hybrid) [56] 


Oil 


Be 


300.7 


300.6(3) [52] 


300.96 (ECG) [53] 


a 2 


Mg 


813.9 


812(6) [52] 




a 2 


Ca 


3,063 


3,081(23) [52] 




Oil 


Sr 


4,645 


4,630(8) [52J 





C. The alkaline-earth atoms 

The use of a fixed core model reduces the calculation of 
the alkaline-earths and their excited spectra to a two elec- 
tron calculation. The two electron wavefunctions were 
expanded in a large basis of two electron configurations 
formed from a single electron basis mostly consisting of 
LTO. Typically the total number of one electron states 



would range from 150 to 200. The use of such large ba- 
sis sets means that the error due to incompleteness of 
the basis is typically very small. The semi-empirical po- 
larization potential needs to include a two-body term to 
deal with the instantaneous interaction between the core 
and the two valence electrons that may be on opposite 
sides of the nucleus [32, [H, H3] ■ 

Details of the calculations used to represent Be, Mg, 
Ca and Sr have been previously described [H, l5lj - [6lf . 
We refer to these semi-empirical models of atomic struc- 
ture as the configuration interaction plus core polariza- 
tion (CICP) model in subsequent text. The matrix ele- 
ment set for Sr incorporated experimental information. 
An experimental value was used for the 5s 2 1 S e -5s5p 1 P° 
matrix element [62^ and the energy differences for the 
low-lying excitations were set to the experimental ener- 
gies 



61] 



III. RESULTS 



A. Polarizabilities 



The static dipole and quadrupole polarizabilities are 
listed in Table |U These are defined using oscillator 
strength sum rules [12, [63[ , namely 



07, 



J On 



(k) 



(2) 



( k) 

where is the oscillator strength for the fc-order mul- 
tipole operator. 

The purpose of Table U is to give an overview of the ex- 
pected accuracy of the present dispersion coefficient cal- 
culations since polarizabilities and dispersion coefficients 
are computed with sum rules using the same oscillator 
strengths. Polarizabilities from high accuracy calcula- 
tions based on relativistic many body perturbation the- 
ory [H, H, H HI, 12, H, |6ll are also shown in Table Q] 
The many body theory results for one electron atoms are 
an all-order relativistic many body perturbation th eory 
with single and double excitations (MBPT-SD) |46l. |65| . 
For two electron atoms, many body perturbation the- 
ory is used to treat the interaction between the core and 
valence electrons while the interactions between the two 
valence electrons are treated with the configuration inter- 
action approach. This is called the CI+MBPT approach 

The agreement with experiment or hybrid experimen- 
tal/theoretical estimates of the static dipole polarizabil- 
ity for Li, Na and K is better than 0.5 %. The level of 
agreement is 1% for rubidium where relativistic effects 
are more important. There is close agreement between 
CICP and MBPT-SD [H quadrupole polarizabilities for 
the alkali atoms. 

There is very good agreement between the CICP and 
CI+MBPT polarizabilities for the light alkaline-earth 
atoms, Be and Mg. Unfortunately, there has not been a 
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TABLE II: The dispersion coefficients (in a.u.) for the ground 
state of alkali atoms interacting with the ground states of 
alkaline-earth atoms. The numbers in the square brackets 
denote powers of ten. Ce coefficients derived from MBPT-SD 
and CI+MBPT dynamic polarizability [§4[ are given in the 
rows with no other C n coefficients and have estimated errors 
given by the numbers in brackets. The notation a [b] means 
a x 10 b . 





c 6 


C s 


Cio 






Be(2s 2 1 S e ) 




Li(2s) 


478.3 
478(3) 


1.743 [4] 


1.038 [6] 


Na(3s) 


522.1 
521(4) 


2. 196 [4] 


1.403 [6] 


K(4s) 


791.4 
790(6) 


4. 775 [4] 


3.924[6] 


Rb(5s) 


869.4 
873(7) 


5.818[4] 


5.076[6] 






Mg(3s 2 1 S e ) 




Li(2s) 


856.8 
853(8) 


5. 676 [4] 


4.535[6] 


Na(3s) 


930.1 

926(9) 


6.727[4] 


5. 735 [6] 


K(4s) 


1,417 
1 411f4) 


1.282[5] 


1.333 [7] 


Rb(5s) 


1,553 
1,556(15) 


1.513[5] 


1.660 [7] 






Ca(4s 2 1 S e ) 




Li(2s) 


1,689 
1,660(14) 


1.417[5] 


1.263 [7] 


Na(3s) 


1,815 
1,782(15) 


1.637[5] 


1.564[7] 


K(4s) 


2,803 
2,756(23) 


3.030[5] 


3.468 [7] 


Rb(5s) 


3,064 
3,030(26) 


3.535[5] 


4.268[7] 






Sr(5s 2 'S 6 ) 




Li(2s) 


2,035 
2,022(3) 


1.898[5] 


1.789 [7] 


Na(3s) 


2,183 
2,167(4) 


2.172[5] 


2.194[7] 


K(4s) 


3,384 
3,362(8) 


3.971[5] 


4.765[7] 


Rb(5s) 


3,699 
3,697(10) 


4.609[5] 


5. 833 [7] 



high precision experimental estimate of the polarizability 
for either of these atoms. However, there has been one 
very accurate calculation of the Be polarizabilities using 
an explicitly correlated gaussian (ECG) basis 53]. The 
ECG polarizabilities should be correct to at least 4 dig- 
its. Agreement with the Be polarizabilities could hardly 



be better. 

Additional information such as magic wavelengths and 
tune-out wavelengths also give information on the polar- 
izabilities of atomic systems [IH, 166l - l72j . Most recently, a 
relativistic variant of the present calculation [7l| was used 
to predict tune-out wavelengths for potassium. Agree- 
ment was achieved with the experimental value (69j . 

The best experimental estimates of the polarizabil- 
ities for Ca and Sr come from hybrid calculations 
where the ns 2 1 S e -nsnp l P° matrix elements from 
photo-association experiments were used to correct a 
CI+MBPT calculation of the polarizability [12|. The 
CICP calculation for Sr, like the hybrid calculations, used 
the experimental matrix element for the resonance tran- 
sition. It is not surprising that it is in close agreement 
with other hybrid theoretical/experimental estimates of 
the Sr polarizability. 

Table U only gives ground state polarizabilities. How- 
ever, excited state polarizabilities u sing the CICP 
method have been given [3(| HH, 0, [58l46l| The polar- 
izabilities of potassium serve as an indicative example. 
CICP calculations gave 615 a.u. for K(4p), 4997 a.u. for 
K(5s), and 1419 a.u. for K(3d) HJ. MBPT-SD calcula- 
tions give 611 a.u. for K(4p 1 / 2 ), 620 a.u. for K(4p 3 / 2 ), 
4961 a.u. for K(5s li2 ), 1420 a.u. for K(3d 3/2 ), and 1412 
a.u. for K(3d 5 / 2 ) UM- The evaluation of the oscillator 
strength sum-rules in previous CICP investigations often 
used calculated energy differences. While this has only a 
small effect on the ground state polarizabilities, the use 
of calculated energy differences will introduce larger dif- 
ferences when applied to excited states where the energy 
differences are smaller and thereby making the polariz- 
abilities more sensitive to small errors in the calculated 
energies. Spin-orbit energy splittings can also impact po- 
larizabilities of the excited states for the heavier atoms. 
For example, the static dipole polarizabilities of the 5pj 
spin-orbit doublet differ by 8% [74j and the present dis- 
persion coefficients should be interpreted as the average 
for a spin-orbit doublet or triplet state. 



B. Ground state dispersion coefficients 

Table HTl lists the Ce, Cg and do dispersion coefficients 
between all combinations consisting of a ground state al- 
kali atom and a ground state alkalinc-carth atom. There 
have been two previous comprehensive tabulations of C n 
coefficients for these combinations of atoms. The tab- 
ulation by Standard and Certain [75j can be regarded 
as obsolete More recently, dynamic polarizabili- 

ties from MBPT-SD and CI+MBPT calculations have 
been used to estimate dispersion coefficients for many 
alkali/alkaline-earth dimers (64[. However these calcula- 
tions were restricted to the lowest order C§ coefficients. 

For all practical purposes, the present C§ parameters 
and the MBPT based Cq coefficients are identical for 
dimers containing Li, Na, K, Rb, Be and Mg. There 
is not a single instance where the two sets of calculations 
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TABLE III: The dispersion coefficients (in a.u.) for the ground state of alkali atoms interacting with the nsnp 1 P° and 3 P° 
excited states of alkaline-earth atoms, n = 2,3,4 and 5 for Be, Mg, Ca and Sr, respectively. The notation a[b] means o x 10 6 . 



Dispersion coefficients 


which could be influenced by an 


accidental degeneracy 


with pseudo-states 


in the alkali atom continuum 


are indicated 


by unde 


r lining. 


















nsnp 1 P° 






nsnp 3 P° 








c 6 


Cs 


c 


c 6 


Cs 


Cio 










Be 








Li(2s) 


S 


1.228 [3] 


2. 370 [4] 


4.737[7] 


4.812[2] 


1.751 [4] 


1.043 [6] 




n 


9.516[2] 


2. 725 [4] 


1.882 [7] 


4.930[2] 


1.785 [4] 


1.061 [6] 


Na(3s) 


E 


1.176[3] 


-3.650[3] 


1.303[7] 


5. 248 [2] 


2. 206 [4] 


1.410[6] 




n 


9.883[2] 


2. 335 [4] 


6.321 [6] 


5.377[2] 


2.250[4] 


1.435 [6] 


K(4s) 


E 


2.173[3] 


7.630[4] 


-5.181[6] 


7.970[2] 


4.807[4] 


3. 948 [6] 




n 


1.611[3] 


8.056[4] 


2.163[6] 


8.161[2] 


4.902[4] 


4.017[6] 


Rb(5s) 


E 


2. 444 [3] 


9. 739 [4] 


5.600[6] 


8.756[2] 


5.861 [4] 


5.109[6] 




n 


1.778 [3] 


1.002[5] 


7. 443 [6] 


8.964[2] 


5.976[4] 


5.199[6] 



Ll(ZS ) 


\ ■ 


Z.oW[o\ 




1 Q AC\ TOl 

l.o4U[oJ 


1 O A A TOl 

l.z44|oJ 


z.oyo|Dj 


o a o i r^i 
2A61[(\ 




11 


1 ooc Tol 

l.ooo[oj 


T 1 A 'Vl A 1 
I .L4:( [4J 


1 A TO \ f 7\ 
lAt Z[(\ 


l.Uoo[oJ 




3. 841) [6] 


AT -/0 \ 


v 

Xj 


Z.Z6\)[Z\ 


4.uyy[oj 


o.yzy[ ( \ 


1.66 ( [o\ 


o q i a rpri 
z.ol4[oj 






r r 
11 


i on^ Tol 

1.o9d[oJ 


—2. 585 [4] 


C ATI 

— 5.971[oJ 


1 1 A o Tol 

1.148[3J 


6.921[4J 


5.053[6] 


K(4sJ 


V 1 


A 1 A o roi 
4. 148 \6\ 


1.300[6] 


i nc7Tol 

i.yD/[8] 


o nf?nrol 
2.0bU[3J 


4.8 (9[oJ 


6.091[<] 




n 


3.131[3] 


2.038[5] 


4.813[6] 


1.760 [3] 


1.443[5] 


1.330 [7] 


Rb(5s) 


E 


4.629[3] 


1.537[6] 


2.606[8] 


2.252[3] 


5.489[5] 


7.382[7] 




n 


3.442[3] 


2.783[5] 


1.870[7] 


1.926[3] 


1.719[5] 


1.708 [7] 


Li(2s) 


E 


-5.557[2] 


1.082 [6] 


Ca 

1.515[8] 


2.307[3] 


6.241 [5] 


6.669[7] 




n 


1.474 [3] 


2.310[5] 


1.714[7] 


2.065 [3] 


1.205[5] 


8. 890 [6] 


Na(3s) 


E 


-2. 476 [3] 


1.281 [6] 


1.742 [8] 


2.438[3] 


6.652[5] 


7.864[7] 




n 


1.096[3] 


3.046[5] 


2. 209 [7] 


2.188[3] 


1.405 [5] 


1.145 [7] 


K(4s) 


E 


1.164[3] 


-9.103[4] 


3.670[8] 


3.871 [3] 


1.151[6] 


1.590[8] 




n 


2. 983 [3] 


-1.946[5] 


6.467[7] 


3.461 [3] 


2. 886 [5] 


2.884[7] 


Rb(5s) 


E 


4.641 [3] 


8.794[5] 


4.611 [8] 


4.223[3] 


1.287[6] 


1.908[8] 




n 


4.099[3] 


1.348[5] 


9. 109 [7] 


3. 778 [3] 


3.442[5] 


3.664[7] 


Li(2s) 


E 


-1.753 [3] 


1.408 [6] 


Sr 

2. 240 [8] 


3.116[3] 


9. 286 [5] 


1.086 [8] 




n 


1.560 [3] 


2.695[5] 


2.235[7] 


2.718[3] 


1.575[5] 


1.257[7] 


Na(3s) 


E 


-5.563[3] 


1.589 [6] 


2.541 [8] 


3.270[3] 


9.838 [5] 


1.266 [8] 




n 


7.231[2] 


3.403[5] 


2.831[7] 


2.860[3] 


1.842 [5] 


1.602[7] 


K(4s) 


E 


3.002[2] 


-2. 382 [7] 


5.054[8] 


5.257[3] 


1.703 [6] 


2.517[8] 




n 


3.428 [3] 


-8.213[7] 


7.551 [7] 


4.582[3] 


3. 769 [5] 


3.962[7] 


Rb(5s) 


E 


1.068 [3] 


4.886[5] 


6.078[8] 


5. 735 [3] 


1.899 [6] 


3.002[8] 




n 


3.927[3] 


-1.103[5] 


9.973[7] 


5. 000 [3] 


4.503[5] 


5. 009 [7] 



differ by more than 1%. There is also better than 1% 
agreement between the CICP and CI+MBPT dispersion 
coefficients for strontium. The largest differences occur 
for the dimers involving calcium, where the CICP C 6 is 
just over 1% larger than the CI+MBPT values which 



use an experimental matrix element for the resonance 
transition. This C§ difference was expected since the 
CICP values of a\ for calcium were just over 1% larger 
than the hybrid CI+MBPT calculation of a x . 
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C. Dispersion coefficients for the alkaline-earth 
excited states 

Table lLUl details the dispersion coefficients for the nsnp 
1,3 po aikaiine-eajth excited states. In the case of calcium 
and strontium, the first singlet excited state is actually 
the ns{n— l)d 1 D e state. For purposes of completeness, 
dispersion coefficients between the ns(n— l)d 1 ' 3 D e states 
states of calcium and strontium with the alkali ground 
states are listed in Table IIVI 

The presence of a downward transition for the 1 P° ex- 
cited state makes it possible for the dispersion coefficients 
to be negative, thereby indicating a repulsive dispersion 
interaction. Contributions of the C n coefficients can be 
negative when the total energy of the transitions originat- 
ing from the two atoms is negative. The 3 P° state does 
not have a spin-allowed transition to the ground state so 
all the dispersion coefficients are positive. Examples of a 
negative Cq coefficient occur for some of the £ states in 
Table HTU 

Tabic Unl exhibits some expected trends. The Cq coef- 
ficients tend to increase as the alkali atoms get larger in 
the Li — > Rb sequence. This is expected since the polariz- 
abilities increase from Li — » Rb. There is also an increase 
in Cq for the nsnp 3 P° states as the atoms increase in 
size from Be — > Sr. This is again a polarizability related 
increase. A steady increase in nsnp 1 P° state Cq values 
does not occur as the atomic size increases from Mg — > 
Sr. These states have downward transitions and the Mg 
3s3p 1 S e polarizability [33] is about twice the size of the 
Sr 5s5p 1 P° polarizability foil ]. 

There are a number of apparent irregularities when ex- 
amining the C n values for a sequence of atoms, e.g from 
Li — > K, or Be — > Sr in Table IIIII and the later tables. 
This occurs because the energy denominators in the sum 
rules for the dispersion coefficients now have single atom 
excitation energies that can be both positive and neg- 
ative. For example, negative Cq values occur for some 
nsnp 1 P° states. What has occurred is that the energy of 
the downward transition of the alkaline-earth exceeds the 
energy increase of some of the upward transitions of the 
alkali atoms. In addition, the energies of the downward 
transition and the upward transition were nearly equal, 
so the energy denominator in the sum-rule was small, 
thereby enhancing the contributions from these terms. 
There were also some negative C% and Cio values for some 
dimers. For the most part, these were also found to be 
caused by a near-zero in the C n sum-rules caused by the 
near cancellation of an energy increasing transition (and 
not necessarily a dipole transition) of the alkali atom and 
the energy decreasing transition of the nsnp 1 P° state. 
In some cases, accidental near-degeneracies in the energy 
denominator leads to dispersion coefficients for one dimer 
that seem to bear little relation to those of another dimer 
for which one would expect similar dispersion coefficients. 
As a specific example, one can refer to the negative C% co- 
efficients for the Ca(4s4p 1 P°)-K(4s) dimer. These were 
caused by the Ca(4s4p 1 P°) de-excitation energy of 0.107 



a.u. being very close to and larger than the K(4s — > hp) 
excitation energy of 0.096 a.u. Something similar occurs 
for the Cg coefficients of the Sr(5s5p 1 P°)-K(5s) dimer. 

There is one potential problem with some of the dis- 
persion coefficients involving excited states when the de- 
excitation energies from the alkaline-earth excited states 
are larger than the ionization energies of the alkali atoms. 
Formally, the energy denominator in the perturbation 
theory sum-rules for these combinations would have a 
zero arising when the excitation energy in the alkali atom 
continuum is equal to the de-excitation energy of the 
nsnp 1 P° excited states. The dispersion interaction in 
this case will have an imaginary part. There is also 
the possibility that an accidental near zero energy be- 
tween the nsnp 1 P° de-excitation energy and the en- 
ergy of the one of the pseudo-states in the alkali atom 
pseudo-continua could lead to an error in the calculation 
of the dispersion coefficients. The only tabulated C n co- 
efficients susceptible to this problem occur in Table IIIII 
The dimers involved contain either the K and Rb atoms 
interacting with the Be(2s2p 1 P°) and Mg(3s3p 1 P°) 
states and the dispersion coefficients that may be suscep- 
tible to this problem are underlined in Table [TTTl 

While the nsnp 1 P° state is the lowest excited state for 
Be and Mg, it is not the lowest energy excited state for 
Ca and Sr. In these cases, the 4s3d 1 D e and 5s4d 1 D e 
states are the lowest energy excited states for Ca and 
Sr respectively. Dispersion coefficients for the 1 D e and 
3 D e states are listed in Table HVl The Cq coefficients for 
the 1 D e states are all positive since these states do not 
have an energy decreasing dipole transition. Some of the 
Cio coefficients are also negative, this occurs because the 
transition energies of the ns — > (n— l)d states are almost 
the same as transition energies of the ns{n— l)d —> ns 2 
transition energies. This leads to the energy denomi- 
nator in the sum-rule used to compute Cio 30] being 
close to zero. The negative C% coefficient for the Na-Sr 
dimer arises from the near-equality of the 3s — > 3p and 
5s4d D e — > 5s 2 1 S e transition energies. 

All the dispersion coefficients involving the alkaline- 
earth ns(n— l)d 3 D e states are positive. This state has 
an energy decreasing dipole transition to the nsnp 3 P° 
state. However, the transition energy for this transition 
is very small and there are no accidental near-equalities 
in energies with any transitions emanating from the alkali 
ground states. 



D. Dispersion coefficients for the alkali excited 
states 

Table fVl gives the dispersion coefficients for the np and 
(?i + l)s alkali states interacting with the ground states 
of the alkaline earth atoms. Table I VII gives the disper- 
sion coefficients between the alkaline-earth ground states 
and the lowest alkali nd states. All the coefficients are 
positive. 

All the C n coefficients in Table [V] are positive and 
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TABLE IV: The dispersion coefficients (in a.u.) for the ground state of alkali atoms interacting with the ns(n— l)d D e and 
3 D e excited states of calcium and strontium. The notation a [b] means a x W b . 





ns(n-l)d 1 D e 






ns{n-l)d 3 D e 




c 6 


Cs 


Cio 


C 6 


Cs 


Cio 



Ca 



Li(2s) 


S 


1.605 [3] 


1.347[5] 


4.943[7] 




3.033[3] 


2.852[5] 


3. 749 [7] 




n 


1.673[3] 


1.356[5] 


1.897[7] 




2. 806 [3] 


1.386[5] 


1.404 [7] 




A 


1.876[3] 


1.357[5] 


1.258[7] 




2.124[3] 


1.203[5] 


1.080 [7] 


Na(3s) 


E 


1.713 [3] 


7.834[4] 


5.944[7] 




2.988[3] 


3.052[5] 


4.215[7] 




n 


1.783 [3] 


1.316[5] 


2.451 [7] 




2. 773 [3] 


1.638[5] 


1.684 [7] 




A 


1.993[3] 


1.567[5] 


1.562[7] 




2.127[3] 


1.392[5] 


1.315[7] 


K(4s) 


E 


2.679[3] 


3.577[5] 


-4.121[8] 




5.330[3] 


6.100[5] 


8.687[7] 




n 


2. 792 [3] 


3.153[5] 


-1.809 [8] 




4.906[3] 


3.328[5] 


3.978[7] 




A 


3. 133 [3] 


3.014[5] 


2.023 [7] 




3.632[3] 


2.861 [5] 


3.095[7] 


Rb(5s) 


E 


2.927[3] 


4.285[5] 


3.549[7] 




5. 840 [3] 


7.103[5] 


1.036[8] 




n 


3.050[3] 


3. 726 [5] 


2.090[7] 




5. 373 [3] 


4.029[5] 


4.920[7] 




A 


3.419[3] 


3. 544 [5] 


4.026[7] 




3. 972 [3] 


3.434[5] 


3.845 [7] 


Li(2s) 


E 


1.996 [3] 


1.228[5] 


8.498 [7] 


Sr 


3.605 [3] 


4.384[5] 


6.361[7] 




n 


2.057[3] 


1.490[5] 


2.591 [7] 




3.312[3] 


1.918[5] 


1.908 [7] 




A 


2.239[3] 


1.658[5] 


1.605 [7] 




2.433[3] 


1.435[5] 


1.275 [7] 


Na(3s) 


E 


2.125 [3] 


-9.792[4] 


9.931 [7] 




3.686[3] 


4.737[5] 


7. 180 [7] 




II 


2.189[3] 


9.393 [4] 


3.324[7] 




3.401 [3] 


2.288[5] 


2. 336 [7] 




A 


2. 380 [3] 


1.908 [5] 


1.983 [7] 




2.545 [3] 


1.698[5] 


1.594 [7] 


K(4s) 


E 


3. 342 [3] 


4.544[5] 


3.521 [8] 




6.271 [3] 


8.915[5] 


1.418[8] 




n 


3.443[3] 


3.853 [5] 


1.453[8] 




5.739[3] 


4.415[5] 


5.377[7] 




A 


3. 746 [3] 


3.669 [5] 


4.833 [7] 




4.144[3] 


3.395[5] 


3. 666 [7] 


Rb(5s) 


E 


3.653[3] 


5.536[5] 


-1.951 [8] 




6.861 [3] 


1.025 [6] 


1.674 [8] 




II 


3. 762 [3] 


4.591 [5] 


-9. 572 [7] 




6.279[3] 


5.286[5] 


6. 642 [7] 




A 


4.090[3] 


4.309[5] 


4. 238 [7] 




4.531 [3] 


4.063[5] 


4.562[7] 



obey predictable trends. The coefficients get bigger as 
the alkaline-earth atoms change from beryllium to stron- 
tium. The coefficients also increase as the alkali atoms 
increase in size from lithium to rubidium. 



The C„ coefficients in Table I VII also obey regular 
trends. One trend is for the Cq coefficients to decrease as 
the alkali atoms increase in size from lithium to rubidium. 
This might seem counterintuitive but the polarizabilitics 
of the lowest alkali nd states tend to decrease in size from 
sodium to rubidium [3(| ■ There is also a trend for the 
C*6 values to increase in size for the larger alkaline-earth 
atoms. 



The Cs and Cio coefficients are negative for the dimers 
with A symmetry. This is not the consequence of transi- 
tions that decrease energy. Rather, the equations |3fjll76j 
for C n coefficients of A symmetry allow the possibility 
that the dispersion coefficients can be negative. 



IV. CONCLUSIONS 

Large scale CI calculations have been used to gener- 
ate polarizabilities and dispersion coefficients for many 
combinations of the alkali and alkaline-earth atoms. The 
underlying accuracy of the calculation is known by ref- 
erence to p revious calculations using the same structure 

model mi m mm. 

The most important coefficients are those for the al- 
kali and alkaline-earth atoms in their respective ground 
states. These coefficients should be accurate to close to 
1% for many combinations and a reasonable upper limit 
of the maximum error would be 2-3%. Somewhat sur- 
prisingly, the only previous tabulation of the full Cq, Cs 
and Cio set for these dimers was the Standard and Cer- 
tain compilation (7o| and the present C n values in many 
cases are an order of magnitude more precise (an assess- 
ment of the accuracy of Standard and Certain compila- 
tion for dimers involving alkaline-earth atoms has already 
been published [Hj]). The same degree of precision is 
not present for dispersion coefficients in dimers contain- 
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TABLE V: The dispersion coefficients (in a.u.) for the excited states of alkali atoms interacting with the ground states of 
alkaline-earth. The notation a[b] means a x 10 b . 









^8 


Oio 














Li(2p) 






T " /O \ 

Li(3sJ 




n.Zii.2 1 oe\ 

tie(As o J 


Zj 


l.ool[oj 


1.00o[0j 


l.D 


4.zzo[oJ 


l.U10[Dj 


z.y / o[oj 




11 


656.4 


1.825[3] 


8.106[4J 








mg(AS o ) 


Zj 


Z.OOU[OJ 


0.04 / [OJ 


/I 7^1 T7l 
4. < 01 [ / J 


7 Qfii r^i 


9 i n i r^i 

Z. 1U1 [OJ 


D.0Zo[oJ 




TT 
11 


i 99/1 r^i 

1 .ZZ4[0J 


A 1 Q7T/11 
4. 10 ( |4J 


9 997Tfil 
Z.ZZ / [DJ 








Col At, 2 1 Q e% l 


Zj 


I . 1 1 1 [Oj 


y.uou[oj 


i i Q/i roi 

1 . 1o4[oJ 


1 . f ou [4J 


4.oUy [DJ 


1 ^c\a roi 
i.ou4 [yj 




TT 
11 


Z.ooU[OJ 


i fiii \k] 

l.Dl 1 [OJ 


o. ( Z0[DJ 








or^os o ) 


Zj 


y. < iu[oj 


1 9^/1 Tfil 
l.Z04[OJ 


i.oy4[oj 


9 9Q1 \Al 
Z. Zo 1 [4J 


D.lo4[Dj 


1 OQ/lTOl 

i.yo4 [yj 




11 


0. ( 01 [OJ 


z.ooy [oj 


1 QS^T7l 
1 .Oo0[ / J 




















Mo (A o\ 

INa^4S j 




i>e(As o ) 


Zj 


Z.Uo4[oJ 


o.ozy[oj 


4.DD0|/J 


4.0 / l[oJ 


1.1 ( o[Dj 


0.0 / y[oj 




n 


1.043 [3J 


a orn Til 

4.350[3J 


2.694|5J 








Mg(6S !3 ) 


Zj 


4.U0Z[0J 


( .40 / [OJ 


i 9nnrsi 

1 .ZUU[oJ 


o.DU0[0J 


Z.4Z4[0J 


7 QQ9TS1 

( .yyz [oj 




TT 
11 


1 G^7TQl 

i.yo < [oj 


0. / OO [4J 


A 91 Q \f<\ 
4.Zlo[Dj 








Hoi A a 2 1 Q e% l 


v 

Zj 


1 1 A SMI 
1 . 14o [4J 


i s i /i r^ii 

l.o!4[0J 


9 OflQ rsl 

z.yuo[oj 


1 099T/11 

i .yzz[4j 


O.Ool [uj 


i.ooo[yj 




TT 
11 


4. f 04 [OJ 


Z.0o0[0J 


1 .00o[ / J 








or^os O ) 


Zj 


1.04 f [4J 


z.ooy [oj 


o Q/iorol 
O.o4o [o\ 


9 QOflT/ll 

z.oyu[4j 


7 nQ9TAl 
( .UoZ[Dj 


z.oo / [yj 




TT 
11 


0.4U0[o| 


4.0OD[0j 


9 A Q 7 rvl 
z.4o / [( J 














T/" /I ™-> \ 

K(4p) 






TY(OS) 




rJe(zs o ) 


Zj 


z.oo4[o| 


o.yoo[oj 


y.oob[/J 


o.obl[oj 


i.yui[bj 


(.out [o\ 




n 


1.413[3] 


9.000[3] 


6.306 [5J 








Mg(dS i3 ) 


Zj 


0.U00[0J 


1 999rfil 

l.zzzlOJ 


9 Q9nrsi 

Z.oZU[oJ 


1. 1U0[4J 


0.o40[DJ 


1 KAA TQl 
1.044[yj 




TT 
11 


Z.0z0[oj 


y.o / o[4j 


fi QOQTfil 

o.oyo[oj 








OolAc 2 1 d e \ 


Zj 


1 9Q1 Ml 
1. Z0 1 [4J 


9 701 Tfil 

z. / yi [oj 


0.010[oJ 


9 A ^QT/ll 
Z.400[4J 


O. / Zl [DJ 


Q coi TqI 

o.ozi [yj 




TT 
11 


0.o4U[0J 


o.zy / [oj 










o„/r „2 1 ne\ 

br^os o j 


Zj 


i.oyo[4j 


o.0Uz[Dj 


7 OQA TOl 

/ .Zo4[oJ 


o.U4o[4J 


1.1Uo[ ^J 


4.4oy[yj 




TT 
11 


/ .oo0[oj 


4.y / o|oj 


O.DD0[/ J 














Kbl^opj 






Kd(ds ) 




n n /i)J 1 ne\ 

rJe(zs o j 


Zj 


z.yoo[oj 


/ .4yo[oj 


1 Ol 1 [Q~\ 

l.oll[oJ 


D.zoU[oJ 


Z.lDo[t>J 


0. /oo[oJ 




n 


1.620[3] 


1.272[4] 


9.225[5] 








Mg(3s 2 'S 6 ) 


E 


5.687[3] 


1.524 [6] 


3.119[8] 


1.181[4] 


4.359[6] 


1.845[9] 




n 


3.009[3] 


1.135 [5] 


8.725[6] 








Ca(4s 2 ^ 


E 


1.363[4] 


3.457[6] 


7.367[8] 


2.624[4] 


9.871 [6] 


4.200[9] 




n 


6.654[3] 


3.797[5] 


2.940[7] 








Sr(5s 2 


E 


1.753[4] 


4.431 [6] 


9.625[8] 


3.256[4] 


1.247[7] 


5.348[9] 




n 


8.331[3] 


5.681 [5] 


4.460[7] 









ing excited states. The excited state calculations involve 
de-exciting transitions that can lead to cancellations in 
the sum-rules used for the computation of the C n coeffi- 
cients. In these cases, the uncertainties in the dispersion 
coefficients can easily exceed 10%. These cases can be 
identified by looking for anomalies in pattern of C n co- 
efficients during an examination over a group of similar 
atoms. 

One limitation of the present calculation is the absence 
of the spin-orbit interaction and the use of LS coupling. 
As mentioned earlier, the difference in the polarizabilitics 



of the Rb(5p) spin-orbit doublet is 8%. This would then 
translate into similar differences in the Cq coefficients in- 
volving this state. Similarly, the spin-orbit energy split- 
ting for the 5s5p 3 P° level of Sr is about 0.001 Hartree. 
Given the close proximity to the 5s4d 3 D e state which 
has a binding energy that is only 0.016 a.u. smaller, 
one could easily see the polarizabilities of the 3 P° spin- 
orbit states differing by 10%, again leading to a differ- 
ence of 10% in the Cq values. The possible impact of 
spin-orbit effects was a primary factor in deciding not to 
extend the calculations of the heavier cesium and barium 
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TABLE VI: The dispersion coefficients (in a.u.) for the lowest nd states of alkali atoms interacting with the ground states of 
the alkaline-earth atoms. The numbers in the square brackets denote powers of ten. 



Ca Cg Cio Cq C's Cio 









Li(3d) 






Na(3d) 




Be(2s 2 ^ 


S 


5.972 [3] 


4.208[6] 


2.273[9] 


5. 849 [3] 


4.073[6] 


2.176[9] 




n 


5.274[3] 


1.408 [6] 


8.463[7] 


5.167[3] 


1.362 [6] 


8.076[7] 




A 


3.181[3] 


-2.547[5] 


-2.674[6] 


3.120[3] 


-2.462[5] 


-2.577[6] 


Mg(3s 2 l S e ) 


E 


1.137[4] 


8.323[6] 


4.581[9] 


1.110[4] 


8.015[6] 


4.389[9] 




n 


1.002 [4] 


2.905[6] 


3.147[8] 


9. 788 [3] 


2.801 [6] 


3.026[8] 




A 


5. 984 [3] 


-2.947[5] 


-1.204 [7] 


5.857[3] 


-2.815[5] 


-1.152 [7] 


Ca(4s 2 1 S e ) 


S 


2.711 [4] 


1.837[7] 


1.033[10] 


2. 564 [4] 


1.747[7] 


9.880[9] 




n 


2.374[4] 


6.338[6] 


8.662[8] 


2. 250 [4] 


6.076[6] 


8.298[8] 




A 


1.365 [4] 


-3.631[5] 


-3.683[7] 


1.307[4] 


-3.435[5] 


-3.379[7] 


Sr(5s 2 ^ 


E 


3.674[4] 


2.416[7] 


1.343[10] 


3.361 [4] 


2.283[7] 


1.282[10] 




n 


3.206[4] 


8.288[6] 


1.221[9] 


2.943[4] 


7.943[6] 


1.168 [9] 




A 


1.805 [4] 


-2.914[5] 


-5.324[7] 


1.690[4] 


-2.835[5] 


-4.573[7] 








K(3d) 






Rb(4d) 




Be(2s 2 1 S e ) 


E 


4.147[3] 


2. 340 [6] 


1.039 [9] 


3.327[3] 


1.646 [6] 


6.464[8] 




n 


3.664[3] 


7.704[5] 


3.625[7] 


2. 943 [3] 


5. 373 [5] 


2.166[7] 




A 


2.215[3] 


-1.396[5] 


-1.272[6] 


1.791[3] 


-9.586[4] 


-8.224[5] 


Mg(3s 2 1 S e ) 


E 


7. 833 [3] 


4.623[6] 


2.125[9] 


6.250[3] 


3.271[6] 


1.337[9] 




l r 
11 


b.yo/ [3] 


l.olofoj 


l.o62[o] 


c Kit? roi 

O.Ol6|3J 


i i a o rd 
1.143[Dj 


i no a rol 
1.034[5] 




A 


4. 129 [3] 


-1.318[5] 


-5.260[6] 


3.315[3] 


-7.370[4] 


-2.896[6] 


Ca(4s 2 ^ 


E 


1.778 [4] 


1.179[7] 


5.300[9] 


1.397[4] 


7.951[6] 


3.094[9] 




n 


1.560 [4] 


4.092[6] 


6.629[8] 


1.227[4] 


2. 782 [6] 


3.263[8] 




A 


9.046[3] 


-9.458[4] 


-1.197[5] 


7.150[3] 


-6.966[3] 


-4.764[6] 


Sr(5s 2 ^ 


E 


2. 226 [4] 


3.871 [7] 


5.877[9] 


1.740 [4] 


1.081 [7] 


4.209[9] 




n 


1.950 [4] 


1.317[7] 


5.303[8] 


1.526[4] 


3.822[6] 


5.827[8] 




A 


1.121[4] 


1.211[4] 


-2. 326 [7] 


8.825[3] 


9.168[4] 


3.038[6] 



atoms. The present results are best regarded as giving a V. ACKNOWLEDGMENTS 

set of average dispersion coefficients. The differences in 

dispersion coefficients involving spin-orbit doublets could 

be magnified if the sum rules contain near degeneracies 
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calculations of the dispersion coefficients to the next level 

of accuracy would require properly relativistic structure This work was supported in part by the Australian 
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